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Abstract 


An introductory review of the large-scale structure and behavior of the 
Earth’s magnetosphere is presented, suitable for inc]usion in courses 
dealing with space physics, plasmas, astrophysics or the Earth’s environ- 
ment, as well as for self-study. The exposition is guided by a series of 9 
quantitative problems, dealing with properties of linear superpositions of 
a dipole and a constant field. Topics covered include open and closed 
models of the magnetosphere, field line motion, the role of magnetic 
merging (’’reconnection"), magnetospheric convection and the origin of the 
magnetopause, polar cusps and high latitude lobes. References and direc- 
tions which may guide readers to more recent developments are also given. 


a 


Prefaoo 


The lost 20 years hove seen oonslderable progress in our understanding 
of the Earth's magnetosphere, the region in which the Earth's magnetic 
field reacts with the solar wind and the interplanetary magnetic field. 

Only occasionally, however, has magnetospheric physics appeared in courses 
on astrophysics, geophysics or plasma physics. 

This article presents on elementary overview of the large-scale 
interactions of the magnetosphere, in a format suitable for inclusion in 
such a course or foi' self-study. To illustrate the coupling between the 
magnetosphere and its surrounding medium, a family of simple fields is used 
— the fields formed by the superposition of a magnetic dipole and a 
constant magnetic field. Three points are given special attention*. 

(1) The subjects considered are of global s^ope, involving the entire 
magnetosphere and its coupling to the solar wind. They prepare the ground 
for any more detailed study of magnetospheric physios and at the sam.e time 
also illustrate basic notions such as magnetic field line motion, merging 
and convection. 

(2) Connections between some of the ideas described and more recent 
research are outlined, and references are given for further exploration of 
various avenues. The added material is mainly for readers who already have 
some familiarity with the subject and much of it may be skipped in an 
introductory study. 

(3) The great virtue of the simple models used here is that many of their 
properties may be derived analytically. Several such calculations are 
included and they may be performed in class or assigned as homework 
(problems 1, 3 and should be solved when they are encountered, since 
tlieir results are used inthe text that follows them). The only 
prerequisites are some basic knowledge of electromagnetic theory and of 
guiding center motion. 
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Since about the year 1000, when the compass was invented, it has been 
known that the Earth has a magnetic field. This field appears to be 
produced by electric currents in the Earth's core and can be viewed to a 
good approximation as a dipole , very nearly aligned with the Earth's 
rotation axis. Since the dipole component falls off with distance at on r 
rate which is slower than that of other internal components, the dipole 
approximation of the internal field improves with increasing distance from 
Earth and will be used in all that follows. 

The geomagnetic field can also be described by magnetic field lines 
C"line.*i of force") which, from the basic definition of magnetic polarity, 
rise from the Earth near the southern magnetic pole and re-enter near the 
northern one (Figure 1). This polarity, incidentally, has undergone many 
reversals during the Earth's geological history, rvQst recently xr> 700,000 
years ago, and it is also opposed to the one observed in the fields of 
Jupiter and Saturn, 

The dipole field may be represented by a scalar potential y 

= - Vy (1) 

and for the geomagnetic dipole, in spherical coordinates (r,0,'f) 

Y = - a B cosO (a/r)^ (2) 

6 

where a is the Earth's radius and B ^30,000 nT (nanotesla, i.e. B is 

equivalent to 0.3 gauss) is the equatorial intensity at r=a. Alternatively, 

( 1 ) 

B^ may also be described in berms of Euler potentials (o»3) 

^ 5 Va X VP (3) 

where 

a = - a B (a/r) sin^6 (^a) 

e 

0 = a? (i)b) 

This representation tends to be difficult to derive unless the field is 


it 


axially symmofcrio, bub ib has bho advanbnge of giving an expllcib 
represenbabion of magtiebic field lines! since 

D.Va = B.vg 5 0 (5a) 

lb follows bhab each such lino is ©harocberizod by poramebric equabions of 
bho form 


a(x,y,z) = eonsb. 0(x,y,z) = consb, (5b) 

At large disbanoes bhe Earth's field becomes negligibly small in 

comparison with bhe i nberplanebary magnebio field (IMF) whidi, because of 

its large scale (compared bo a), will be assumed bo bo a constant field 

The IMF is of solar origin and is carried away from bhe Sun by bhe action 

solar wind (SW), a steady radial outflow of hob plasma from bho 

Sun's corona. At the Earth's orbit a typical velocity for the solar wind is 

-3 

V s 300 km/s and a typical density is ns 10 cm ! for convenience, a local 
corbesian frame will be used in which the x axis points sunward and the 
northward half of the dipole axis will be assumed to coincide with bhe 
positive z axis. Actually, bhe angle between bhe dipole and bhe sunward 
direction may be as much as 35° out of bhe perpendicular, because bhe 
Earth's rotation axis is not perpendicular to the ecliptic and is also 
inclined to one dipole axis, but this complication will be neglected here. 

A typical magnitude for the interplanetary field ^ is around 
5 nT, and due bo the nature of bhe interaction between the SW and the IMF 
(see further below) while B tends to be smaller than either of 

the other components. In what follows, however, ^ will not be constrained. 

As an approximation of fields which bend to ^ as r-^O and to as r+®, 
we shall study simple superpositions of the two fields 


B 


o 


( 6 ) 


It turns out that such combined fields divide space quite naturally into 
2 regions — an inner region dominated by which will be named 
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(2) 

(following Gold ) the magnetonphere , and an outer region dominated by D^. 
Many important features are not reproduced by such crude models, including 
the collision-free bow shook created in the SW at a distance <r2a upstream 
of the magnetosphere; the values of v and n cited above are indeed more 
appropriate to the region behind the shock, while ahead of it v tends to be 
larger and n smaller. 

Furthermore, contributions to B from electric currents flowing in the 
magnetosphere are ignored, including those from the long current sheet 
stretching the magnetospherio "tail" on the night side (Figure 2), and many 
plasma processes are grossly simplified. Nevertheless (as will be seen) 
many basic phenomena are iiiusbrated quite well at this level. 

The Closed Model 


Consider first the possibility of a purely northward IMF 


~o o~ 


(7) 


^ may be represented by a scalar potential 


-<0 -- - 


( 8 ) 


or also by axisymmetrical Euler potentials 


a = (aB^/2) (r/a)^ sin^O 


(9) 


$ = ai 


( 10 ) 


The scalar potential for the combined field is obtained by superposition 

Y = Yo + Yd 01) 

For Euler potentials, in general, superposition does not work, but if (as 
here) two fields share the same a superposition of a gives a correct 
result. Thus in the present model 
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B B Va X VB (12) 

a s a sin^o C(BQ/2)C»Va)^~ B^Ca/r)) (13a) 

B B aY (13b) 

PROBLEM 1 ; 

(a) Show that the surface a=0 is a sphere of radius 

I /% 

r„ B a (2B-/0_) (Where else in space does a vanish?). 

0 6 0 

By (5a), B is everywhere tangential to this sphere, hence 

it completely encloses the dipole (Figure 3). Show by means 

of (11) that B„=0 on r = r_. 

r o 

(b) Show by means of (11) that B vanishes only at two points, 
located at the Intersections between the sphere and the z 
axis. Points at which B=0 are termed neutral points . 

This model will be called "dipole in sphere"; the sphere natu- 
rally forms the boundary of the magnetosphere or magnetopause 
and, at least on the side facing the Sun, resembles the boundary 
actually observed. The observed distance to the magnetopause 
10a) is somewhat smaller than r^ derived here, because the 
pressure of the solar wind plasma, which has been ignored, pushes 
the boundary further earthwards. This generally also produces a 
finite current flow along the boundary. 

The "dipole in sphere" is the arohtype of closed models of the 

magnetosphere which assume a completely enclosing magnetopause 

through which no field lines pass. It bears some resemblance to 

Hill's spherical vortex in ideal fluid dynamics and to a two- 

(3) 

dimensional model of plasma filaments studied by Schindler • Tn 
the magnetospheric context it was discussed at least as early as 

(H ) 

i959 I and a more realistic "closed" model was advanced in 1960 
(5) 

by F , Johnson , who has presented arguments in favor of a 


c* 
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variant of the "olooed” magnetosphere as reoently os 1978^^ 

The "ciipolo in sphere" oan be modified to aeoomodato any 
direction of B., If for instance B_ shifts to the direction of 

•Hy *==V 

the unit vector bi it is only necessary to rotate the external 
field pattern of Figure 3 until its asymptotic direction lies 
along b. Or in other words, one then defines a new system of 
spherical coordinates with the polar axis along Jb and then assu- 
mes that the potential y oV (11) 

2 

Y s -a COSO (a/r) “ cosO (It) 

is given in the old system for r<r^ and in the new one for 
r>r^. At rsr^ the field is discontinuous, indicating a current 
flow along the magnetopause, but as was already noted, such 
currents indeed occur in nature. An example of such a closed 
model, with b = X, is given in Figure t. 

Closed models will be further discussed below, but at this 
point we shift our attention to the outer region with r>r^ and in 
particular to the question, how can the existence of the IMF be 
reconciled with the flow of the SW. 
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Tho Elecfcflo Field 

Far from bho origin both the raagnotlo field and the SW velocity have their 
undisturbed interplanetary values 

B vs- v^l (15) 

Suppose first that at any point all charged partiolos of the SW have the same 
velocity. How then can they make any headway in a direction orthogonal to D, 
rather than be forced into tight little spirals around the z axis? Tho answer 
is that an electric field E is set up 

E = - V X B (16) 

so that bho magnetic force q(^xB) on any particle with charge q is completely 

cancelled by the electric force qE, allowing the particle to proceed 

undisturbed along the x axis (for details about how E is set up, see for 

( 7 ) ~ 

insbsncs chapt. 4 of the text by Longmire ), More precisely, eq. 16, also 
called the "magneto-hydrodynamic (MHD) condition" must be modified to include 
non-electromagneblc forces (e.g, gravity) and terms representing 
accelerations, but all added berms are generally negligible under the 
conditions encountered here, save perhaps in regions very close to a neutral 
point (see further below). The MHD condition only gives the component of E 
orthogonal to B; it is customary to add the requirement 

E„ = 0 (17) 

since charges generally can easily slide along magnetic field lines and thus 
equalize any potential difference which violates (17). Thus (16) indeed gives 
the total electric field, and by (15) 


E = - 

(18a) 

$ s yv„D^ 
•'0 0 

(18b) 


If SW particle velocities cover a finite range, the two above forces in 



gynorni do nofc balaneo, The sifcuabion may then be viewed in bho frame of bho 

SW, in whioh bho (non-relabivtsbic) local eloobric field vanlohoo by (16); in 

bhob frame mosb parbicleo have a amall rosiduoX veXociby and bhus carry oub 

some small-scale gyrabion, bub aparb froin bhls bhey share bho voloeiby v^ of 

bho frame. Or else, viewed from bho frame of bhe Earbh, guiding eonbor bhoory 
( 7 ) 

requires bhab all porbiclos (regardless of charge) share an eloobric drlfb 

C19) 

which sabtsftes (16), From either poinb of view, the bulk voloeiby of the 
plosma remains bhe same. 

This eonsbanb recbillnear velocity is modified, however, in bho vioiniby of 
rsT^. Inbuibively, one would expeob bhe flow bo divide up in the manner of 
Figure '1, where now bhe lines subside rsr^ are regarded as strortinlinea rather 
bhan as magnetic field lines. To see whether this indeed happens, one proceeds 
as follows. 

The electric potential 4' can be viewed as a function of (x,y,z), or 
equivalently of 3 independent functions of (x,y,z): let two such functions be 
(a,p) of (13). For bhe third one could take s, the distance along a field line 
to the point in question, measured (say) fran the plane z=0: the exact choice 
is immaterial, in any case, for by (17) 

B,V4> s 0 (20) 

It then follows from (5a) that no matter what the third coordinate is, * 
does not depend on it but on (a,P) alone. This of course is consistent with ♦ 
being constant along a field line, as implied by (17). 

In order to find the electric potential near r=r^ (and from it, through eq, 
19, the bulk velocity) it is thus sufficient to evaluate the function ^(aiP) 
far away from the origin and then propagate it along field lines to regions 
closer to the dipole. At large distances bhe second term of (13a) may be 
neglected, leaving 
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a n {Bq/2o)(xV> s (B^/2a) y^/sln^CB/a) C8D 
From (18b) then 

«(o,P) 5 Vq C2aBg)^^^ 0 ^^^ sln(B/a) (2?) 

Tf (2B) is valid at any one point on a field lino, it is valid anywhere else 
on it, and it will therefore also give ^ near rsrg, provided the full 
expression (13a) is used for a. It is now straightforward to derive E and then 
use (19) to find the flow velocity at any point in the field. If only flow 
lines are required, however, a faster procedure exists, since (19) implies 


v.Vy = 0 

(23) 

v.y* s 0 

(2‘t) 


Hence streemXines are characterized by lines along which y and 4> ore constant. 


PROBLEM 2 ; Derive the equations 

(a) Of a streamline in the (x,y) plane which crosses the y axis 
at a distance (5/*0r^ frcxn the origin, 

(b) Of a streamline In the (x,z) plane which crosses the z axis 
at a distance (5/*0r from the origin. In either case the 
app.-oximation 5 <i' 2 may be used. For each line determine: 

(1) Does the shape resemble one of the contours of Figure 

(2) Does the shape match one of those contours? 

(c) What flow lines, if any, follow the surface of the sphere 
rsr^? If such lines exist, what will be the flow velocity 
along them? 



Oonvflctton 
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What flow voloeifcy exists inside the closed magnetosphere? The present model 
offers no oluo •»- indeed, the internol flow is eomplotely deeouplcd from the 
external one and o porfeetly aeeeptable solution for r<r^ would bo 

V 9 e E a 0 (35) 

It is easy to visualize, however, that if any plaame fills the roeion r<r^, 
and if maoentum is tronsferred to it oeross the boundary, then just inside the 
boundary the plasmo will tend to flow in the (-x) direetion (’’antioMnward"), 
However, steady fluid flow must be oontinuous, so that if a steody-stabe 
situation is established, a return flow in the (•♦■x) direction must exist 
Domowhero in the interior { o possible pattern for such a flow, as viewed in 
the (x,y) plane, is shown in Figure 5a, This motion is often eallea 
maKnatospherio oonveotion , because it resembles, qualitatively., the convective 
flow of a heated fluid (Figure 5b). 

Such a flow, driven by a "viscous like" momentum transfer across the 
boundary, was indeed suggested in the earliest theory of magnetosphsrio 
convection, proposed by Axford and Hines^° . Their reasoning proceeded in the 
opposite direotion from the, one given earlier; they interpreted the motion of 
auroral ares and related phenomena and concluded that magnetospheric 
disturbance-s were often associated with arrival of plasma from the night side, 
around the midnight meridian, This, they suggested, was the return flow from 
and antisunward plasma motion near the boundary, and they proceeded to unify 
both flows in a global theory of magnetospheric convection. 

The concept of magnetospheric convection remains one of the fundamental 
ideas in all theories of the magnetosphere. However, few things in nature are 
simple: as noted below (and indeed, as already recognized by Axford and Hines) 
there also exists an alternative explanation of convection, based on the 
so-called open magnetosphere . 

The Dunaeysphere 


If the field in (7) is reversed to represent a "purely southward" IMF 


a difforcnfc axtoynmotrical model of the inogneboaphorc io obtained. Ifc ia 
called hero the ’'Dungoyophero", bocauae it was firab proposed by James Durgey 
in 00 on archbypo of bhe "open'' moBnoboaphorio model, further 

described below. A cross section of the Dungeyophere io ahown in Figure 6 and 
its propertloa are readily derived by an extension of bhe calculation 
performed ear' ier on the "dipole in sphere" model, The dobailo are liobod in 
Problem 3 and the derivation io loft os on oxcoroiee bo the reader, though 
final results will be oiled. 

PROBLEM 3 t For a "Dungoyophore" produced by the superposibion 
of bhe oipole field of (D-O) and the consbonb 
southward field of (36) 

(a) Derive the scalar potential y. 

(b^ Derive the Euler potentials (a, 8), given asaCr.O), Usaf , 

( 0 ) As Figure 6 suggests, bhe equatorial plane of the Dungey- 
sphore contains a circular "neutral lino" along which DsO. 

Find the radius r^ of this line and the value which a 
assumes on it, and compare r^ to the radius of the spherical 
magnetosphere produced with a constant northward field 
having bhe same intensity 

(d) Using «Q, find the colatitude 0^ on the surface of the 
Earth from which magnetic field lines lead bo the neutral 
line. It is permissible to neglect B^ in comparison with B^. 

If = 5 nT and B^ a 30,000, what is bhe value of 9^? 

(e) Derive bhe electric potential ^Ca,8), using bhe same method 
as in eq. (32). 

(f) What is bhe electric field E in bhe Earth's ionosphere 
(assume r°a), in (x,y,z) coordinates? With B^ = 5 nT, v^ = 


300 Icm/s, whot is chs largest voltage drop between two points 
in the polar area bounded by 

Answers: 


(a) Y = a oosO CB^Cr/a) - B^Ca/r)^] (2?) 

(b) a = -a sin^o [(B^/2)(r/a)^ + Bg(a/r)3 (28) 

(c) Eq.(2r) has the form y = a cos6 f(r), and on the neutral 
line all derivatives of y must vanish. The factor ooso 
assures the vanishing of 3Y/3r anywhere in the equatorial 
plane, but the other derivative only vanishes if f(rQ)=0. 

This gives 

r^ s a (29a) 

1 /3 

a distance smaller by a factor 2 than the spherical 
radius in the preceding model. From (28) 

= -(3/2) aB^ (B /B^)^‘^^ (29b) 

V \J V w 

(d) Neglecting the external field, at r = a 

Oq = -aBg sin^Gg (30) 

Substitution of (29b) gives 

0Q « 16.7° (31) 

(e) The derivation follows that of (22), except that the sign of 
a is reversed, so that a is replaced by (-a), a positive num- 
ber, Also, because ^ is reversed, the sign of ^ in (18b) is 
reversed, giving 

♦ (a, 3) = -Vq (23Bq)^^^ sin(3/a) (32) 
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(f) Neglecting B. in comparison with B_» at rsa 

Q U 

-o = aBjj sin^O 

« = -Vq (2BpBg)^''^ y (33) 

The polar electric field is thus a constant field, directed 
from dawn to dusk. The largest voltage exists between extreme 
points on the dawn-dusk meridian, where OsO^ and sin 'f = ±1, 

Its value in volts is found from (33) if MKS units are used, 
and amounts to 

= 2v^a (BB^Bg)^'^^ sin0g ^ 6 10^ volt (3^) 

The most significant qualitative property of the Dungeysphere is that it is 
"open"; while in the "dipole in sphere" only two types of field lines exist, 
purely terrestrial and purely inter planetary, the Dungeysphere admits a 3rd 
type, so-called open field lines, connected to Earth at one end and extending 
into space at the other. All open field lines connect to circular patches with 
0<0g (or O>n-0g) around the magnetic poles, which will be termed the "polar 
caps." In these caps a constant electric field exists, directed from dawn to 

5 

dusk and with a predicted voltage drop of about 6 10 volts. 

Surprisingly, most features of this crude model agree with observations. 

Roughly circular patches of approximately constant dawn-to-dusk electric field 

have been observed near the magnetic poles by various spacecraft and 

( 15 ) 

also by other methods (see sect. 8 of a review on the subject ). The 
boundaries of the patches roughly agree with the predictions of the model and 
they may be monitored by observing the locations of quiet auroral arcs 

( 1 6 17 ) 

("auroral oval") which tend to parallel them a few degrees equatorward ’ 

An electric field across the polar cap implies an electric current across 
the polar ionosphere, completed by field-aligned current sheets, flowing 
downwards into the caps' dawn boundaries and upwards from the dusk boundaries. 
Such sheets, carrying 1-2 million amperes, have indeed been 
observed^^®*^^’ they are flanked by weaker sheets of opposite polarity, 
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located further equatorward and apparently having a diffferent origin. 

These currents may he visualized as arising from a "classicai” dynamo, from 
a closed electrical circuit in which (whatever frame one uses) part of the 
current flows in a stationary medium and part in a moving medium which crosses 
magnetic field lines. Viewed in the frame of the Earth (see ref. 15, Fig, 18) 
this circuit would consist of the ionosphere (stationary), the solar wind 
(moving) and, linking the two, magnetic field lines anchored at the dawn and 
dusk edges of the polar caps. 


The most obvious discrepancy between the model and observations is in the 
magnitude of A4>, which falls short of the derived value by a factor of 10 or 
more. Perhaps the reason is that the actual magnetosphere, due to the pressure 
of the solar wind, is greatly elongated in the x direction, while the 
Dungeysphere is axisymmetric. Thus in the Dungeysphere the polar sheaf of open 
field lines has a circular cross section and is proportional to its asymp- 
totic width Ay, in the y direction; in the actual magnetosphere, this sheaf is 

model's value 

elongated in the x direction to perhaps 10 times the^ , and if its 

magnetic flux remains unchanged, its width in the y direction, and 
consequently A^, must shrink by a similar factor The above may sound as a 
plausible explanation, bub it is still remarkable that in spite of the great 
deformation of the outgoing sheaf, the size of the polar cap and the direction 
of E there differ so little from the prediction of the Dungeysphere. So far, 
however, no more detailed theory of open field lines exists. 


PROBLEM 4 : If E = what is the shape of the interface 

between the Dungeysphere and the solar wind? 

Answer : By equations (23)- (24) streamlines of an ideal MHD flow are 
distinguished by the constancy of y and of ♦, though here only the former is 
required. By (27) and (29a), y vanishes on the neutral line, and therefore the 
surface bounding all solar wind streamlines will be the surface ysO, which is 
a sphere of radius r^. Thus, figure 6 notwithstanding, the Dungeysphere, too, 
is a sphere I 


PROBLEM 5 : Let an axisymmetric surface current be added along the 
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surface rsr^j (see preceding problem) such that the added 
field produced by it is 

for r<r^ 2 = const. (35) 

for r>r^ s -7^1 = dipole field 

= (r^jB^/2)(r^/r)^cos 0 (36) 

(a) If the combined field B is represented by toroidal and poloidal "poten- 
tials” '1'^ and vi»g, through 

B = Vx'l»^£ + ^xVx'l'211 

find ijfp (i|i. vanishes) . 

Hint: 

2 

VxVxtligll 5 V(8/Dr)\li2r - £ V \|i2 

In source-free regions the second term vanishes, allowing tjig to be related 
to the scalar potential y- 

(b) Derive a(r,0) for the combined field, for regions Inside and outside rsr^. 
Hint: ^ may be expressed in two ways 

A 

B s Vx(axp) = Vx o(a/rsino) f 

B^ = Vx(Vx(})2£) = 7x(V'!»2x£) = -7xC(atli2/30) f 3 

(c) (optional) Using a computer-controlled graphical device, draw field 
lines of this model in the (x,z) planes (electric potentials 

and other properties may also be derived). This simulates, 
at least near the noon meridian, the extra current needed for 
achieving pressure balance at the boundary, when the pressure 
due to the solar wind plasma is also taken into account. 
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Mogneblc Merging and Convesfclon 

The curl of the MHD condition 

E 5 -vxB (37) 

gives 

3B/8t - VxC^xB) = 0 (3B) 

/ p n p 1 \ 

This may be shown ' '(the proof is given in many plaana texts and is 
completely analogous to vorticity preservation in ideal fluid dynamics) bo be 
equivalent bo two conditions; (1) "flux preservation", the requirement that 
the magnetic flux threading a closed line of particles, each moving with the 
local bulk velocity, remains unchanged in time; and (2) "line preservation, 
that two such particles which initially share a magnetic field line, continue 
to do so at all times. The second condition (which is weaker than the first, 
and may be deduced from it) tends to be especially useful in visualizing MHD 
flow, and it readily leads to the notion of "moving field lines." Suppose a 
given field line is defined at b=b^ by the particles strung out along it; 
then the field line defined by the same particles at any other time may be 
regarded as one and the same field line, having "moved with the plasna." 

An application of this is provided by problem 2-c. Consider a field line 
distant enough from the dipole so that it may be considered bo be parallel bo 
the z axis, and let it cut the x-axis at (x,y) s (Xq, 0). Particles located on 
this line at b=b^ and having a sufficiently large jz] will pass through one 
of the branches of the z axis at some later time tg, which means that for an 
arbitrarily small At, at tg-Ab they are sunward of the axis and at tg+At 
anti sunward of it. 

But what about particles with small jz] which shared the same field line at 
b=t^ and which are therefore destined bo hit the sphere r=r^ at t=bg? At 
tg-At they are upstream of the sphere, at tg+At they are downstream of it, so 
that in an arbitrarily short time 2At they must have skirted the sphere with 
an arbitrarily large velocity (the side of the sphere around which the 
particle chooses bo move is only defined if its initial y at t=t^ is 
displaced from zero an infinitesimal amount to one side or the other). Though 
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this eonftlusion moy appear surprising it is borne out by explicit 

-1/P 

caloulotioni since £ contains 34>/Da which by (22) is proportional to a t 
and a-*-0 as n-^r^. 

As a more significant example, consider 5 solar wind particles, emitted 
with a radial velocity from the Sun's equatorial corona one day opart and 
numbered in the order in which they are anitted. At time tsO, just before the 
first particle emerges, they are all clumped together (Figure 7a) and thus 
certainly share the same magnetic field line, Five days later they are all 
moving radially in space, and presumably are still joined by the some 
magnetic field line, but that line is a spiral rather than a radial one, 
because the Sun rotates about 15° between successive emissions (Figure 7b), 
This spiral shape (with large fluctuations) has been confirmed by 
observations and makes an angle of about ^15° with the radial at the Earth's 
distance, which is the reason for the earlier assertion that in the IMF v' 



Nevertheless, "line preservation" runs into trouble with the open 
magnetosphere (Figure 6). Consider two particles, numbered "1" and "2", 
sharing the same field line as it approache'/ Earth. A short time later "1" is 
on an open field line above the northern c»p while "2" is on an open line in 
the opposite hemisphere, and the two no longer share the same field line. On 
the other hand, each now shares its line with a collection of ionospheric 
particles, though no such sharing existed before the particles encountered 
the magnetosphere, or is likely to exist after they have passed beyond its 
influence. 

The breakdown occurs at the moment when the field line containing the 
particles passes the neutral point N1, and indeed theory shows that the 
concept of line preservation must be modified here, that it does not assure 
the continuation of "sharing a field line" beyond points where the line 
passes a point with BsO. The process involved here is called "magnetic 
merging" or "reconnection" and in it the interplanetary field line AB flowing 
into N1 is cut into two parts. Simultaneously, a "closed" field line CD 
flowing into N1 (in a manner described further below) is also severed, and 
each half of the terrestrial line links up bo half an interplanetary line to 
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form on "open” Xine, Later, ob a neutral point N2, the process is reversed •— 
the two interplanetary segments ore reconnected and continue their outward 
journey, while the terrestrial sections ore also reunitec and flow bock 
earthwards. 

(2?1 

Hiis idea of ''magnetic merging" has led bo many theoretical studies , a 

few laboratory experiments and some attempts to observe it in space, but a 
more complete discussion is well beyond the scope of this introduction. 
Although Dungey^^^^ and others have regarded merging os a process in 
which magnetic energy is given up accelerate and heat porticlei (especially 
in solar flares), observational evidence suggests that if merging does occur 
at the front of the magnetosphere, it does not involve the release of very 
much energy. As for the nighbside reconnection at N2, observations are 
scarce: the consensus is that if this takes place, its location is well 
beyond the orbit of the moon ('/'60a) and that it probably does not proceed in 
a smooth and steady fashion. 

The flow of closed field lines into N1 snd out of N2 suggests that 
convection is unavoidable in an open magnetosphere, unlike what is found for 
a closed configuration, where (25) is a valid solution. This indeed is the 
case, as can be seen by considering (2H) , which requires drifting particles 
to follow surfaces of constant potential From among such particles let us 
give particular attention to an ionospheric ion which is held close to Earth 
by gravity and which keeps a roughly constant altitude h of (say) 200 km. By 
(33) the projection of a polar equipotenbial surface on the surface h=200 km 
is a line of constant y, extending between the two points at which it cuts 
the polar cap boundary 0=0^ (northern cap). Presumably, the ion follows this 
line (Figure 8). However, equipotential surfaces cannot end suddenly: a 
fringing pattern must exist at 0>0 q, i.e. on closed magnetic field lines, in 
a way suggested by the broken lines in the figure. 

This fringing electric field 1| will cause particles such as the one we have 
selected to continue their drift 

V = ExB/B^ 
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beyond bhe boundary 0=0 qi uud beaause of (2iJ) i the lines whieh such parbioles 
follow will preaarve their values of Tlvus the ionospheric portiole at 
position ”1" in Figure 8 will bend to move along its equipobenbial surface to 
position "g" , then (on an open magnetic field line) bo "3"» out of the polar 
cap bo "<1" and finally back bo "1". However, because of line preservation, 
all particles which share a field line at ”1" will move together to ''2“ , a 
position corresponding (very nearly) bo line CD in Figure 6. Ttien after 
merging, since line preservation still applies to all particles that pass on 
bhe same side of BsO, half of them will stay together and will ultimately 
share a line with "3". while the other holf will move above bhe opposing 
polar cop. Finally, at "Jl" both halves ore united again bo form a line 
similar bo EF of Figure 6, after which they conveob once more sunward. 

All this is qualitative. In particular, the above argument does not predict 
whether bhe fringing pattern is pressed close to the polar cap boundary or 
extends far from it. In the first ease, bhe sunward return flow in bhe 
equatorial plane is pressed close to bhe magnetopause and does nob penetrate 
close to Earth, while in bhe second case bhe penetration is much more 
pronounced. A complicating factor nob mentioned so for is the existence near 
Earth of an additional electric field E^, caused by the Earth's rotation and 
tending to impose co-rotation upon the plasma: the interested reader is 
referred bo a more complete review for details. 

A useful mathematical model for the fringing pabbbern is obtained as 
follows. By (32), the polar electric Held (0<0^) has the form 

^■(a.e) = “-tg (a/a^)^^^ sin(e/a) (39a) 

(since a^<0, the negative sign may be omitted). A class of fringing fields 
which fit continuously with bhe above is then 

4>(a,s) = (oQ/a)*^ sin(p/a) ( 35 b) 


PROBLEM 6 : 


(a) Compare the magnitudes of bhe two berms in (28) at r=rQ, 
near the equatorial plane., Which is the larger one, and by 
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what ratio? How does this ratio change os one advances in- 
wards from 

(b) Show that if the smaller term in 08) is neglected and if 

(39b) is valid with Ksl, a constant dawn-to-dusk electric field will 
exist on closed magnetic field lines in the equatorial plane. 


While theory at this level does not predict the fringinj, field, it may be 
deduced from observations, either close to the lonosphere^^’^^^^^^ or near the 
equatorial plane . All such analyses seem to support (39b), at least os a 
rough app''oximation , but suggest that the value of K is not 1 but 2. 



Asymmetrle Models 

Deeouse the DunBeyspliGre is nxisymmetrie, B in it vonishes along an entire 
line. However, in more general situations, where sueh symmetry no longer 
exists, B usually vanishes only at two points '(or at a larger even 
number). For instance, if the field configuration is stretched on th^ night 
side (as observations suggest), only two neutral points remain, analogous to 
N1 and N2 in Figure 9: this development was anticipated in the discussion of 
the merging process, which thus remains unchanged. The detailed 3~dimensional 
geometry of magnetic field lines in the vicinity of an isolated neutral point 
Is not always easily describable on a qualitative level' ' • 

A different example of an open magnetospheric model with ? isolated neutral 
points, studied in problem 7, arises when the dipole field of equations 
(l)-(2) is superposed upon an arbitrary IMF 




PROBLEM 7 5 Derive the coordinates of the neutral points formed 
when the constant field ^ of C»I0) is superimposed 
upon the dipole field of equations (1)-(2). 

Note: The problem may be studied in cartesian coordinates 
(x,y,z) and if B >f0, appropriate intermediate variables are 
Gsx/y, Usz/y and r"^. At a certain stage a quadratic equation 
arises and one might be led to believe that each of its solu- 
tions represents one of the neutral points, but it isn't so. 

(a) Show that one of the above solutions is non-physical, while 
the other gives both neutral points, when all choices in 
assigning algebraic signs are utilized. 

(b) Solve for (x,y,z) given BQy=0« 



Uiiless («0) represents a purely northward field (os in eq. 7) the resulting 
model is olwoys open* The "polor eop" stays roughly cireulnr even when when 
Bq„bO, but its size shrinks steadily oa the direebion of approaches that 
of the z axis, all of which agrees in a quallbobive way with 
observations . There even exists a slight sunward sliift of the center of 
the polar cap in response bo whieh also is in qualitative oeeord with 
observations The structure of the electric field predicted by these 

models, on the other hand, is quite intricate, and its details probably do 
nob fib observations boo well, 

( O 1 N 

The preceding results must all be evaluated numerienily by computer » 
because (a,0), and hence ♦(o,6), can no longer be derived onalytioally. 
Problem 8 below contains hints for deriving ♦ for such models: the task is 
considerably lengthier than the solutions of other problems listed here, but 
is still well within the range of a graduate projeot, Program listings in 
FORTRAN (IP sec. on a VAX 11/780) or BASIC (5 hours on an HP ?830) are avail- 
able upon request. 

In the actual magnetosphere, lb should be noted, bounriary currents always 
exist and with such currents, both open and closed magnetosphere models may 
be constructed. If the models ore open they may resemble those described 
above, or else they may have boundary currents such as those appearing in 
problem 5- If they are closed they always have boundary currents (unless 
is strictly northwards) and resemble the configuration of Figure with the 
external curves representing magnetic field lines and having an arbitrarily 
directed asymptotic direction, 

PROBLEM 8 : Given a magnetic field whieh is the superposition 

of (il0) and a dipole, map the electric potential in 
the polar cap. 

Procedure: 

(a) Define a rectangular grid of 21x21 (or more) initial refe- 
rence points, centered at the magnetic pole. A convenient 
approximation is to place the initial points not on the cur- 
ved surface of the Earth bub on the plane z= tangential bo 



ifc ot tho magnofcle pole, A spoeing between the points equi- 
volenb to 1-2 degrees oi' latitude usually ossures that the 
entire polar eop is Itioluded and an orroy PCI.iJ) whieh will 
store (ot the end of the ooloulotion) the oleetrio potentials 
of the points should bo defined. 

(b) Define a reference plane orthogonal to of (JI0), ot a 
distance Rq“30q to ROa from the origin. On that plane find 
tho point closest to the origin, ond assign to it 

on orbitrory reference potential (e.g. zero). Use (UO) to 
derive the electric potential ot large distonees and apply 
the result to express tho potential of on orbitrory (x,y,z) 
on the reference plane. (Cout'on! there exist two reference 
pianos, one threaded by field linos from the northern polar 
cap ond one linked bo the southern cap. Be sure the one 
selected is the one oppropriabe bo tho cap being studied!) 


(o) Taking in burn each of the reference points of port Co), 
brace numerically the magneti',., field line which rises from 
it, using a sbepsize v'O.Sa (initial steps may be smaller , 
later ones larger). Allow a certain number of steps ant -neck 
each time if the reference plane has been crossed. If a cros- 
sing has oecured, derive the crossing point (linear inter- 
polation is sufficiently accurate) and calculate its poten- 
tial; by (17), this is also the potential P(I,J) of the point 
at which the tracing began. If the field line fails to reach 
the reference plane in the allotted number of steps or inter- 
sects Earth, diagnostic outputs such as 0 or -1 should be 
assigned to P(I,J). 


(d) Print out the array P(I,J) in suitable units with 2-3 figure 
accuracy and interpolate manually bo get the polar boundary 
and equipobential contours (or else, program the computer to 
perform the task). One interesting study is the variation in 



tho size of the polor cop ao a function of the direction of 
D^, allowing tho magnitudo to atay eonstant^^^\ 


Polar Cusp and Toil Lobes 

There exists another type of asymmetry which Is of great practioal interest 
in the actual magnetosphere- As noted in problem ili if all SW partieles hove 
the same velocity and obey the MHD condition, then tho boundary of tho 
Dungoysphere is spherical. In actual fact, however, SW particles always hove 
a finite range of velocities and of directions of motion, and this spread is 
greatly increased by their passage through the bow shock of the Earth. Thus 
the spherical boundary (redrawn in Figure 9) will not impose on absolute 
barrier to their motions instead, as soon as merging has occured, some SW 
particles will begin spilling over it onto the terrestrial port of newly 
merged "open" lines. 

Because particles flowing into a stronger magnotic field are subject to 
magnetic mirroring, most such particles will not penetrate very for and will 
quickly move out again, all the time staying on the same moving field line as 
other solar particles which had shared their line earlier bub hod not moved 
earthwards, A few select porbioles among those that penetrate, with veloci- 
ties almost parallel to B, may however roach as far down as the ionosphere. 

By the above argument one might expect a region of .solar wind plasma 
reaching all the way to the ionosphere and shaped, in profile, like the 
curved funnel cloud of a tornado. Such a region is in fact observed^^^^"^^^^ 
and it is known as the polar cusp . 

If the external field is purely southward, as in (26), both polar cusps are 
completely equivalent, More frequently, however, the IMF exhibits the spiral 
structure discussed earlier (in connection with the motion of magnetic field 
lines), so that while the SW velocity ^ far from Earth is radial, its 
components v„ and Vj. , parallel and orthogonal to B, are approximately equal. 


Now to maintain a SW motion orthogonal to B (i.e. , to maintain Vx ), the 
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clcefcrle field of (37) Is weeded i ond ib is this field • too, that defleefca 
bhe flow around the tnognefco sphere. On the other hand, beeauso of (17) t v„ is 
not affcebcd (bo lowest order) by the eleebrio foroe. This suggests a eerbaln 
asymmebry between the polar caps, beeause depending on bhe mognebie polarity 
of interplanetary field lines ("away" from bhe Sun or "toward'' the Sun) open 
field linos from one of bhe oops will extend in the direobioR of bhe Sun 
while those from the other oap will point away from it. If v„ moves SW 
parbielos along bheir guiding field lines like beads along a wire, one would 
expeeb plasma bo bo shoved deeper into the "sun-faolng" cusp ond thus eroabe 
an asymmetry between the two polar regions, 

Many observed asymmetries eorrelabe with B^y of ('JO), which is generally 
taken as an indicator of the magnebie polarity of the In 

particular, when energetic portieles are created by solar flares near oebive 
sunspot groups, they often hove higher intensity (at first) at bhe sunfacing 
polar cap than in opposite polar region^^^^^**^^ . Not much is known as yet, 
however, about unequal penetration of plasma related to Inberplanebary field 
polarity, though some striking correlabioris between bhe intensity of cusp 
electrons and hove been reported ^ . 

Beyond the cusp, in the regions of open field lines extending to the 
downstream end of bhe magnetosphere (ore Q-N2 in Figure 9) the situation is 
quite different, These regions are called bh® "high latitude lobes" (or "tall 
lobes") and in bhe actual magnetosphere they are quite large, since the 
difference in x between Q and N2 is stretched bo a length about 10 times 
larger (or more) than in t‘ s Dungey sphere model. In this region the solar 
wind is moving alongside the magnetosphere or (mostly) away from it, and only 
very few of Its particles manage to buck the general trend and flow 
earthwards, On the other hand, terrestrial particles can easily escape along 
open field lines in these regions. 

One would therefore expect the plasma density n in the lobe region to be 
remarkably low, and this is indeed observed. In the solar wind niA>10 cm in 

the cusp n»/’2-3 cm , in the magnetosphere 1 cm is typical (though close bo 
Earth n becomes much larger) and in the "plaatia sheet" in the middle of the 
tail, near zsO, n,/’0.3 cm"^. In the lobes, however, nvrO.01 cm”^, making this 


c# 
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by far the moat rarefied plasma in the Earth's neighborhood. 
Assessment 


Which model oomes closest to reality? Is the magneto “sphere open or closed? 
The preponderance of the evidence so far points towards an open 
magnetosphere; the polar electric field t the convective flow, the relation of 
the size of the polar cap to B„_ (and one may add, the strong correlation 

(Up) 

between Bqj, and magneto spheric disturbances ), all these suggest that the 
magnetosphere is open. So do the existence of polar cusps and tail lobes, and 
phenomena which correlate with B^y. 

The situation is not always clear-cut. If the magnetosphere is closed and 
convection is driven by a viscous-like momentum transfer, the polar pattern 
of electric equip- tentials will still resemble Figure 8, at least 
qualitatively. However, it will not be easy to explain why the polar caps 
have rather abrupt boundaries, or why their size correlates with B^^* 

Again, the existence of ci sps is not conclusive evidence, because it was 
pointed out^^^^ that in a closed configuration (like that of Figure 3) the 
neutral points would tend to cave in to the pressure of the external plasma 
and this would also produce cusp-like features. 

Even if an open magnetosphere is accepted problems remain. In a dipole 
field, and also in the Dungeysphere, the lowest magnetic intensity B on any 
field line is always in the equatorial plane. This property is important in 
the theory of trapped radiation, where the minimum of B on a field line acts 
as a natural midpoint, across which trapped particles bounce back and forth. 
For the dipole- in- sphere, however, this holds only up to a certain distance. 
Closed field lines which cross z=0 not boo far from r=r^ will have their 
smallest values of B near their points of closest approach to the neutral 
points, and at z =0 they will actually have local maxima of B. 

PROBLEM 9 : Given that a field line crossing the equatorial 

plane of the "dipole in sphere" in the range 
Xro<r<ro has there a local maximum of B, what is > ? 
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Hint; Derive B . In the equatorial plane, aB /ao=0, showing that B has 

an extremum there. At rsXr^ the extremum switches from minimum to 

maximum; what can then be said about the second derivative? 

This property of a minimum turning into a local maximum is common to many 
analytical models proposed as approximations to the actual magnetosphere, 

(uiiwiJKS 

from the Chapman-Ferraro model' to more elaborate image dipole 

representations'**^^ and models utilizing external harmonic s^ **"^^ • 

Theoretically, it has been noted^**^^ that the longitudinal ("second”) 

adiabatic invariant may vary discontinuously where this happens. What seems 

surprising is that studies of average field intensities^^®^ and of pitch 

(51) 

angle distributions of trapped particles in the Earth's magnetosphere 
suggest that the actual magnetosphere also has this property, in contrast 
with the Dungeysphere which does not. 

A different problem involves the electric field pattern of Figure 8. While 

observations have tended to confirm it, there exist indications of a 

different structure during prolonged spells of northward with not 2 but 

4 "convection cells" in the polar caps^^^^^^^\ The meaning of this is not at 

all clear and it is hoped that the twin "Dynamics Explorer" mission, slated 

(54) 

for launch in 1981 , will gather enough data on this puzzling effect for 

appropriate magneto spheric models to be constructed. 

Do there exist any alternative models to those described here? The closed 
model has 2 types of field lines, the Dungeysphere 3, and a model with 4 
types was obtained by Podgorny et al.'^^'^ in laboratory experiments which 
simulated plasma flow around the magnetosphere. These workers found that each 
"interplanetary" field line and each "terrestrial" one, upon merging, split 
not into 2 parts as in the Dungeysphere but into 3; the northern and southern 
sections then joined up as in a simple open model, while the two middle 
sections formed an independent eddy or "visor" (Figure 10) , or possibly a 
chain of such eddies^^^^. A variant of this model may well come closer to 
. existing conditions than any of the simpler models, but observations have not 
yet produced conclusive evidence on this point. 
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Conoluslon 

We have attetnpbed fco present an elementary introduction to magnetospheric 
dynamics, including open and closed configurations, flow patterns, electric 
fields, field line motion, convection, merging, cusp regions and tall lobes, 
using some simple analytical models. The results may be compared bo Figure 2 
which represents a schematic view of the actual magnetosphere. Apart from the 
largo day-night distortion the simple models have provided a fair qualibBtive 
fit to observations, though it should be realized that the physios of the 
plasma sheet, the substorm and aurora, the radiation belt and its "ring 
current" are all beyond such crude approximations and require additional 
study, 
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CAPTIONS TO FIGURES 

Figure 1 The geomagnetic dipoie. 

Figure 2 — Schematic view of the actual magnetosphere (Mizera 
and Fennell, Rev . Geophys . Space . Phys . , 1<47, 1 978) 

Figure 3 — The "dipole in sphere" produced by superposition of 
a constant northward field and a dipole. 

Figure R — (no caption) 

Figure 5a — Convective flow in the equatorial plane of a model 
magnetosphere, schematically shown. 

Figure 5b — Convective flow in a heated fluid. 

Figure 6 — The " Dung ey sphere" formed by superposition of a con~ 
stant southward field and a dipole. The particles 
marked "1" and "2" illustrate the violation of line 
preservation caused by magnetic merging. 

Figure 7a — Five particles initially located in one region of 
the solar corona, illustrating how the property of 
line preservation can provide the shape of inter- 
planetary field lines. 

Figure 7b — The 5 particles of Figure 7a and their shared mag- 
netic field line, 5 days after the first among tliem 
was emitted from the Sun. 

Figure 8 — Electric equipotential lines in a map of the nor- 

thern polar ionosphere, in an idealized schematic 
view. The numbers give successive positions of a 
particle (and therefore also of the footprint of a 
magnetic field line) participating in a convective 
flow. 

Figure 9 — Illustration of the polar cusp and of the high 
latitude lobes, using the Dungeysphere. 

Figure 10 — Podgorny's magneto spheric model, with R kinds of 
field lines. 
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APPENDIX; S 0 1 u b i o n s 


S e 1 e 0 b e d 


Problems 


b 0 


Problem 1 ; « „ 

(a) a s a sin‘=^0 C(B^j/2)(r/a)‘^ - Bg(a/r)] = 0 

gives (r^/a)^ s SB^/B^ or ® 

By (2), (8) p 

Y = - a eosO CB (r/a) + B (a/r)*^] 

B vanishes if « 

3y/ 3I' = “ COSO [B - 2B.(a/r)'^] s 0 

O 6 

This again gives bhe oondibion for r^.. 

(b) As shovfn above, B^ vanishes on bhe sphere o=0. The obher componenb Bq 
is proporbional bo 

3y/ 30 = a sinO [B_(r/a) + B„(a/r)^3 

O 6 

For r>0 bhe brackebs are always posibivo, so bhis berm vanishes only if 
sin0=0. 


Problem 2 ; 

Tal A sbreamllne confined bo bhe (x,y) plane has aubomabically y=0 and 

bherefore (24) has bo be used: 

$ = 0 

and by (22) „ 

o sin'^f = 0 

2 

At bhe crossing poinb (x,y) = (0, 1,25g ), sin f =i, and since bhab poinb is 
also in bhe equaborial plane where sin 9=1, by ( 13 a) 

a = a[(B /2)(a/r)C (r/a)^ - (2B/B ) ] = a 

o e o o 

3 

Subsbibubing r = (5/4 )r^ and using bhe approximabion (5/4) <^2 shows bhab bhe 
firsb berm in bhe brackebs is aboub bwice bhe size of bhe second one. Thus 

aBg (a/r) = aB^ (4/5)^'^^ 

Since bhe flow line is complebely confined bo bhe equaborial plane sin0=1, 
ibs equation may be wribben „ _ 

[(B^/2)(r/a)‘^ - B^(a/r)3 sin I’ = 

which indeedghas bhe same form as (13a), burned around by 90 ° so bhab sin^y 
replaces sin‘^0 . 

(b) A sbreamline confined bo bhe (x,z) plane has aubomabically 4>=0 and 
bherefore (23) should be used; 

2 

Y = -a COS0 [B^(r/a) + Bg(a/r) 3 = Yq = consbanb . 



Yq = “ a (o/r)^ C(r/a)^ + -SoBg (a/r)^ 

Yq - (l6a/5) 

The flow line's equation is thus approximately 


COSO [Bjj(r/a) + Bg(a/r)^3 s (16/5) (B^Bg/»0^^^ 

It does not match the field line equation (13a) in any way. In fact, at 
large distances it reduces to 

z/a = (8/5)(r^/a) > (5/4)(r^/a) 

Hence the flow line decreases its distance from the x axis as it approaches 
the origin, in contrast with the external field lines of Figure <1. 


Problem 5 ? ? 

(a) Let B = -Vy t Y « a cosO [g°(a/r)‘^ + g°(r/a)3 


by 


(all this is easily generalized for an arbitrary harmonic function). Then 
Y = - (8/3r)(<|>2r) 


dTg will have a corresponding expansion, involving the same harmonic functions. 
Equating coefficients « o 

ilig = a oosO tg°(a/r)‘^ - (l/2)g°(r/a)3 

Noting that the external field is southward , one obtains: 

For r<r^ g° = -B^ 8° = 

For r>rp g° = -B^ + (r^j/a)^(B^/2) g° = B^ 

(b) 

o = - sino (r/a) (3<li/80) 


From this (again, this is easily generalized for any axisymmetrical 
a = a sin^o [g°(a/r) - (1 /2)g°(r/a)^3 

One then substitutes the expresions derived in part (a). The above provides 
an alternative derivation of (13a) and (28), or equivalently, of (9). 

Problem 6 ; 

(a) Equation (38) is p p 

a = - a sin 0 C(B^/2)(r/a) •+■ B^(a/r)3 

If the brackets are opened, this is resolved into 2 terras a.+a_, and their 
ratio is 

r 

“/“2 = (B^/2Bg)(r/a)^ 

3 

In the Dungeysphere (r^/a) = B^/B^ , so that at r=r^ the ratio is 1/2 and 

the second term is twice the first one. As one moves inward the disparity only 



A-3 


increases, because diminishes while Og Bfows. 

a 

If one approximates a - a sin r B (a/r) 

V 

p 

then at sin Osl, oq. (39b) gives 

« s - (♦ /B )(-« )( r/a) sinl* = - («_/B )(-a )(y/a) 

0 C v> QUO 

Since (““q) is positive, this gives a constant dawn-to-dusk field, in the 
direction of vy. 

Problem 7 ? 

Ti) Let all distances be measured in units of a. The dipole field may 

be written « n 

B^j VYjj with Yj = “ " '^e 

If the constant external field has components (B^.Bg.B,), the cartesian 
components of the total field have all to vanish at the nedtral points 



= ^1 - 

3BgXz/r^ = 0 

(1) 

"y 

= Bg- 

3B^yz/r^ = 0 

(2) 

«z 

= B3- 

- 2 2 ‘i 

(3z - r^^)B /r-^ = 0 
6 

(3) 


The first two 

eqs. may be rewritten 




1/r^ = (B^/3Bg)/xz 

(4a) 



1/r^ = (Bp/3B^)/yz 

(4b) 

assuming B„/0 

and combining 




or 

x/y = 

(5a) 



X = Gy 

(5b) 

This allows x 

to be eliminated from 

(3) 



Bg - (Bg/ByzXBz^- x^ - y^) = B^ - y^(1+G^)]/3yz = 0 (6) 


Let a new variable U and a new constant H be defined: 


U 5 z/y 

(7) 

H . B3/B2 

(8) 


3H - 2U + 

(1+G^)/U = 0 

(9) 


- (3H/2)U - 

(1+G^)/2 = 0 

(10) 

U = (3/4 )H + 

[ (9/4)H^ + 2(1+G^) ] 

1/2 

(11) 


Two solutions exist, of opposite signs. The negative solution however is 
non-physical , because it implies opposite signs for y and z and therefore , by 
(4b), a negative value of r^. Hence only a single solution exists, with U>0. 
Multiplying (4b) by r 
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1/r^ = (Bg/3B^)(x^+ y^+ z^)/yz s (B^/3B^J)C(1^<}^)/U + U] (12) 

which allows r bo be derived. Other variables follow readily, e.g. for x 

r^ = x^+ y^+ z^ s x^Cl + (1+U^)/G^3 (13) 

Note that if (x,y,z) is a solution, (-x,~y,-z) is one too, because it loads to 
identical values of (G,U,r), the intermediate variables form which everything 
follows , 

If BpsO, (2) gives ysO at the neutral points, but (^la) still holds. 
Substituting in (3) 

B- - (B,/3xz)(2z^ - x^) = 0 (Iti) 

Let now ^ 

U' = z/x (15) H' = Bo/B, (16) 

Then „ ^ 

(U')‘^ - (3H'/2)U> -1/2 = 0 (17) 

U' = 3H'/»1 ± [9(H')^/M + 2]^^^ (18) 


Again, only the positive solution is acceptable, and from it 

1/r^ = (B^/3Bg)(U» + 1/U') (19) 

Two symmetric solutions result from this, one with (x,z) and the other with 
(-x,-z) . 

Problem 9 *. 

Y r. - a COS0 CBg(a/r)^ + B^(r/a)] 

B = COS 0 [2B (a/r)^ - B_3 = cosO F(r) 

r 6 0 

B 0 = sin 0 [Bg(a/r)^ + B^] s sin 0 G(r) 

„2 2 „ p 2 2 „ p 2 

B s cos 0 F + sin 0 G 

3 B ^/30 = 2 coso sine (G^ - F^) and vanishes if 0=it/2, cos 0 =O. 

= 2 (cos^r - sin^e)(G^ - F^) 

If 0=n/2, this vanishes only when G = ±F. If G=F then 

B (a/r)^ = 2B 
e o 

r/a = (B /2B^)^''^ = (r /a) 

e o o 

which yields the required condition. If G=-F we get (a/r)=0 and hence no 
solution exists for any finite r. 
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Pi’oblom 8 ; Llnting of the matn program, In FORTRAN. 


i*» » AS* > BD I F » SI 


CttMMON^ XO . YO I ZO . R. , R7 , BX I » BY I , BZ I , B I , A I T A 
XiiYl.Zl 

DIMENSION P(m*>,SO)>FC:2Z*> 

1001 FORMAT <SF©.3) 

1 002 FORMAT ( ©F0 . 0 ) 

1003^F0RMAT^(^j BX.PY.BZ OF IMF«‘* r SFO, Si OAMMAi DIPDLE« ’ , FS. Oi 

1004' FORMAT ( ^ STEP IN RE. •% F8.3> BIST. TO REF FL=-%Fi?..0/) 
lOOS . FORMAT < ^ .^POLAR STEP« : , FS, 3, •' DEG., |W SPEED £ ’ , Ni. 6, 


1006 

1007 


::o 


c 


1008 


40 

50 

100 

110 

150 

200 


FORMAT < . ' central^ VALUES OF ( X , Y ) «*' , 2F© . 3/ ) . 

READ (5,1U01) BXI»BY1»BZI ! coffipononts of imf 

RtAD (oilUOl) '.sLi’pP 1 SL in RE is stCfP alons field lino 

negative into north pole. SP is polar 
step in decrees. 

READ (5,1001) XIOiYIO ! Central values of plotted part 

of polar cap. To plot entire cap, 
u 0 luG 2or*ci *Por* both 

READ (5,1002) BDIP, VSW, RP, VOLT ! dipole term in samma, SW 

speed in km/s, unit in which volt- 
ages are measured (2000 v OK) 

WRITE (6,1003) BXI,BYI,BZI, BDIP 
WRITE (6,1004) SL, H9 
WRITE (6,1005) SP, VSW, VOLT 
WRITE (6,1007) XIO,YIO 

BI c, SORT(BXI#BXI BYI«-BYI + BZI#BZ1) ! IMF intensitv 
RP ~ (-1. )#F!9 ! for north mas pole dist is nesative 

DO 20 iJ‘-*l,22 
FjJ) a J-1 

1,22) '.labels tabulation 
' (X9,Y9,Z9) are coordinates of 

i point a distance RP from the 
! orisin with r alons B 


uk 


WRITE (6,1006) (F(J), 

X9 K R9 # BXI/ BI 
Y9 « R9 -IS BYI/ BI i 

Z9 = R9 * B2I/ BI ! 

PI - 3.1415926 
•SP » SP <«• PI/ 180 ! 

SS >=6.371 ! 

VO >= VSW S5/ VOLT ! 

VO (BZl 

DO 30 J1 = 1, 21 
DO SO J2 = 1, 21 
P(wU,J2) - 0. 

DO 200 II « 1, 31 
DO 150 12 = 1, 21 
XO => XIO + SP (11-15) ! 
YO >= YIO + SP # (12-10) ! 
ZO « 1. 


polar step in radians 
earth radius in 1000s of 
IMF potential = 

^ (Y-Y9) - BYI * (Z-29)) 


km 


(XO,YO,ZO) is initial point 
in the polar cap 


XI2 « XO 
XLO = BXI 




13 


L2 « 1 
DO 100 
SI = SL 
XI3 = FLOAT 
IF (XI3.lt. 21) 
IF (XI3.LT.6) 


(Xa-X9) +• BYI # (Y0-Y9) + BZI (Z0-Z9) 

XLO is criterion for crossins the reference 
Plane, where it chanses sisn. 

! _^L2=1, no crossins vet, L2=2, plane crossed 

if 

I 

( 13) 

! Initial steps are half sir*? 
SI ~ SL-K-0.25 ! ..or quarter sice. 


200 

Re* 


'.standardize step size 


Radial dist. , f9. 6) 


add 100 


CALL ADVANCE 
IF (R7.GT.1.0) GO TO 40 

^format ( •'^ ( 1 1 , 12, 13)«-', 315, 

P(I1, 12) = -2, 

GO TO 150 

_ CALL REFPLANE ( XLO, L2, X9, Y9, Z9) 

GO TO (100, 50 ) , L2 

P(I1,I2) = (VO* (BZI*YO - BYI*ZO)) -f- 100. 

PI ane . crossed, evaluate Potential, 

^ ^ to eliminate nesative values. 

GO TO 150 
CONTINUE 

P(I1,I2) = -1. ! Sisnifies ref plane not reached 

steps. 

CONTINUE 

WRITE (6,1006) XI2, (P(I1,K), K=1 , 21) 

TYPE 1006, XI2, (P(I1,K), K=l, 21) 

CONTINUE 

STOP 

END 


in 200 



Problem 8 ; Listing of subroutines, In FORTlVkN 


SUBROUTINE BFIELD 


•derives, unit vo'Ztor (A1 jA2iA 3) 
alons the direeticn of the tielid 


COMMON XO, YO, ZO> R2. R7> BXI , BYI i B7 1 » Bl , A1 , A2, A3. BDIP, SI 
1, Xl.Yl.Zl 

R2 XO#XO + YOifrYO + ZCmO 
R7 s S0RT(R2) 

RS c R2tfR2*R7 

B1 ^ BXI 4* 3 .« BDIP ZO XO/RS 
B2 « BYI + S.* BDIP * ZO * YO/RS 
B3 a BZI 4*C3.'»> ZO ZO ~ R2)# BDIP/R5 
BO a S0RT(B1*B1 4- B2-frB2 h- B3K-B3) 

A1 a Bl/ BO 
A2 a B2/ BO 
A3 a 03/ BO 
RETURN 
END 


SUBROUTINE ADVANCE 


advances a step alons field line 


COMMON XOt YOt ZO, R2. R7,. BXI , BYI , BZI , BI , Al , A2. A3. BDIP, SI 
Xl,Yl,Z.l ' 

CALL BFIELD 

XlaXO 

Yl=YO 

Zl«ZO ! save i;<1d point 

XO a XO 4- S1#A1 


ZO a ZO 4* Sli«-A3 ! first approx., new point. 

IF (Zl.ECi.l.) CO TO 10 ! at 1st step no old pt. exis^ 
CALL BFIELD 

XO « 0.5 (XO 4 - XI 4 - SH^A1) 

YO a 0.5 (YO 4 * Y1 4 - SH«-A2) 

ZO = 0.5 » (ZO + Z1 + S1U-A3) 

final approx to new point 

RETURN 

END 

SUBROUTINE REFPLANE ( XLO. L2, XS, YS. Z9) 

COMMON X0,Y0,Z0,R2,R7,BXI,BYI,BZI,BI, Al.A2,A3,BDIPrSl 
Xl.Yl.Zl 

XLi a XLO .'save old value 

XLO a BXI*(X0-“X9) 4 - BYI*(Y0“Y9) + BZHM Z0-Z9) 

XM a XLO XLI ! negative when XLO changes sign 

IF (XM) 1,2,3 

S2 a Si ! save oU value 

SI a SI* ABS( XLI/ (XLI - XLO)) ! last step to brim 


! at 1st step no old pt. exists 


XO a XI 
YO a yi 
ZO a Zl 

CALL ADVANCE 


' XLO)) ! last step to bring 
us (almost) to reference plane 


XLO a 0. 
SI a S2 
L2 a 2 
RETURN 
END 


Restore standard step 



